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ABSTRACT 

We report the first result of weak gravitational lensing survey on a 2.1 deg 2 7? t -band image taken with a wide 
field camera (Suprime-Cam) on the prime focus of 8.2 m Subaru Telescope. The weak lensing mass reconstruction 
is applied to the data to search for dark matter halos of cluster scale; M > 10 14 M Q . The reconstructed convergence 
field is divided by l-cr noise to obtain the signal-to-noise ratio map (S/N-map) of the detection. Local maxima 
and minima are searched on the S/N-map and the probability distribution function (PDF) of the peaks are created 
to compare with model predictions. We found excess over noise PDF created from the randomized realization on 
both positive and negative sides. Negative peaks imply the presence of voids in the dark matter distribution and 
this is the first report of the detection. Positive peaks, on the other hand, represent the dark matter halos and the 
number count of the halos on the 2.1 deg 2 image is 4.9 ± 2.3 for S/N > 5 where the Gaussian smoothing radius 
of the convergence map is 1'. The result is consistent with the prediction assuming the Press-Schechter mass 
function and the NFW halo profile under the cluster normalized CDM cosmology. Although the present statistics 
is not enough due to the limited field of view, this work demonstrates that dark matter halo count via weak lensing 
is a promising way to test the paradigm of structure formation and cosmological model. 

Subject headings: cosmology: observations — dark matter — large scale structure of universe — gravitational 
lensing 



1. INTRODUCTION 

Dark matter halos, which are the virialized systems formed 
via the gravitational amplification of initial density fluctuations, 
are one of the most important cosmological probes. Since their 
evolution is governed purely by gravity, their number density 
and the evolution of clustering properties can be precisely pre- 
dicted analytically as well as numerically. Weak lensing is a 
unique observational technique to detect dark matter halos be- 
cause it detects halos not by their light but by more essential 
physical parameter, the mass. It also allows us to directly es- 
timate the halo mass bypassing a use of ambiguous empirical 
relations like the mass-temperature (or velocity dispersion) or 
mass-luminosity relations. Therefore, comparison of theories 
and weak lensing observations is a direct and reliable way to 
test the paradigm of the structure formation and to constrain 
the cosmological parameters. 

Theoretical prescriptions for predicting statistics of weak 
lensing halos have been developed so far. Kruse & Schneider 
(1999) first calculated the expected number of dark matter ha- 
los detected by means of the aperture mass method (Schneider 
1996). Predicted number density of halos exceeds 10 per deg 2 
for cluster-normalized cosmologies if one obtains moderately 
deep images whose galaxy density is n g ^30 arcmin" 2 . Bartel- 
mann et al. (2001) pointed out that the number density is sen- 
sitive to the density profile of halos. On the other hand, Jain 
& van Waerbeke (2000) discussed the statistics of the peaks on 
convergence map. They claimed that the statistics can be an 



alternative approach to probe the mass function of dark matter 
halos and to constrain f2 m . 

Observationally, weak lensing has become reliable measures 
to map the dark matter distribution in clusters of galaxies (see 
a review by Mellier 1999). Recently, serendipitous detections 
of mass concentrations outside clusters on a reconstructed mass 
map have been reported by several authors (Erben et al. 2000; 
Umetsu & Futamase 2000; Wittman et al. 2001; Miralles et al. 
2002). All cases have no clear optical counterparts except the 
one by Wittman et al. (2001) which is identified as a cluster of 
galaxies at z = 0.276. 

These findings demonstrate that dark matter halo search via 
weak lensing is feasible and can detect even lower luminos- 
ity dark matter halos that are overlooked in conventional opti- 
cal cluster surveys. The number of samples obtained so far is, 
however, limited. In order to make statistically significant dis- 
cussions, weak lensing survey over a representative cosmic vol- 
ume is crucial. The survey requires not only wide field but also 
high image quality to make reliable lensing analysis. Suprime- 
Cam on the 8 m Subaru telescope is the most ideal camera for 
the purpose. The median seeing in I c band monitored over a 
year and a half is reported to be ~ 0.6 arcsec (Miyazaki et 
al. 2002). In this Letter, we report the first results from the 
Suprime-Cam weak lensing survey that covers contiguous 2.1 
deg 2 field. We examine statistical properties of local peaks on 
the reconstructed convergence map by comparing with model 
predictions. 
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2. OBSERVATION AND DATA REDUCTION 

We set field size at 2.1 deg 2 , which is the largest size pos- 
sible to be accomplished during the commissioning phase of 
the Suprime-Cam. The field we chose is centered at R.A. = 
16*04 m 43 s , decl. = +43° 12' 19" (J2000.0). In order to compare 
light and mass we observed one of the seven fields that Gunn 
et al. (1986) made cluster surveys (GHO survey). The 2.1 deg 2 
field has twelve GHO clusters and one Abell cluster (A 2158). 
Note, however, that the number density of clusters in this field 
is comparable to the one averaged over the seven GHO fields. 

We obtained /? e -band images on the nights of 2001 April 23- 
25. Suprime-Cam provides a field of view of 34' x 27' with a 
scale 0".202 pixel" 1 . Nine contiguous fields were observed in 
3x3 mosaic configuration. Each exposure on a given field was 
offset by 1 ' from other exposures to remove cosmic rays and 
defects on the CCDs. Total exposure time on each field was 
1800 sec (360 sec x 5). 

The individual images were de-biased and then flattened us- 
ing a median of all object frames taken during the observing 
run. Registration on each field was first performed assuming 
a simple geometrical model. The parameters of the model in- 
clude (1) optical distortion, (2) displacement and rotation of 
each CCD from the fiducial position and (3) the offset of point- 
ing between the exposures (Miyazaki et al. 2002). These pa- 
rameters can be solved by minimizing position errors of control 
stars common to all exposures and used for the first transforma- 
tion of the images to a common coordinate system. After the 
transformation, the residual of the position of the control stars is 
<~ 0.5 pixel rms. The residual is due to several effects not con- 
sidered in the simple model including atmospheric dispersion 
and asymmetric aberration of optics. We note, however, that 
the residual vector is smoothly changing over position and can 
be well modeled as the third order bi-linear polynomial func- 
tion of position. This model then gives the fine correction of 
the transformation function. As a result the final residual de- 
creased typically down to 0.07 pixel rms (14 miliarcsec). The 
seeing in the resulting image is 0".68 FWHM and the scatter 
among the fields is quite small of 0".04 rms. 

We used the software suite imcat, an actual implementation 
of Kaiser et al. (1995), to carry out object detection, photom- 
etry and shape measurements of objects. Catalogs created for 
the nine fields are registered using stars on overlapping regions 
to result in a final catalog whose total field of view is 1 .64° x 
1.28°. Differences in photometric zero point among the fields 
due to variation of sky condition were compensated at this stage 
but the adjustment was not significant (<~ 0.05 mag). For weak 
lensing analysis, we adopt galaxies of 22 < R < 25.5 whose 
signal-to-noise (S/N) ratio, nu, calculated in the imcat exceed 
15 (Bacon et al. 2001). The number of objects in the catalog is 
297547 (39.3 arcmin -2 ) and the count peaks at R = 25.2. 

3. GALAXY SHAPE ANALYSIS 

The estimate of the shear, 7, from the observed ellipticities of 
galaxies, e = -hi^hi} /(hi+hi) involves two steps where 
Iij are quadrupole moments of the surface brightness of objects. 
First, PSF anisotropy is corrected using the image of stars as 
references, 
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scattered over the field of view and we made fifth order bi- 
polynomial fit to the values as a function of the positions. This 
function is used in Eqn.( 1) to correct the ellipticities of faint 
galaxies. The nine fields are treated independently in this pro- 
cedure. The rms value of ellipticities of the reference stars, 
(|e*| 2 )2, is reduced from 2.7% to 1.3% as a result of the cor- 
rection. Note that the rms before the correction is already small 
thanks to the superb image quality of Subaru telescope. 

Luppino & Kaiser (1997) worked out how to correct the el- 
lipticities for the effect of the seeing. The pre-seeing shear 7 is 
described as 



7 = p:V', p y = p sh - 
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where P sm is the smear polarisability tensors and is mostly di- 
agonal (Kaiser et al. 1995). (P* m )~' e * was calculated for stars 



where P s h is the shear polarisability tensor. The P 7 of individ- 
ual galaxies are, however, known to be noisy estimate and we 
thus adopted smoothing and weighting method developed by 
van Waerbeke et al.(2000; see also Erben et al. 2001 for a de- 
tail study of the smoothing scheme). For each object, twenty 
neighbors are first identified in the r g -magnitude plane where 
r g is a measure of size of objects which the imcat yields. A me- 
dian value of P 1 among the neighbors is adopted for a smoothed 
P 7 of the object. The variance of raw 7 before the smoothing 
among the neighbors, cr 2 , is used to estimate the weight of the 
object, w n , as w„ = l/(er 2 +a 2 ) where a 2 is the variance of all 
the objects in the catalog and a ~ 0.4 here. Under the weight- 
ing scheme, an averaged value of a certain observable, (Q), is 

calculated as J2f=i Wn.iQi/ J2f=i w n,i instead of J^f =l Qi/N. The 
(tr(P 7 )) for all objects is calculated as 0.31 but P 7 depends 
on r g . For objects with r g > 3.5, (tr(P 7 )) is almost constant, 
<~ 0.65, and becomes small below the r g down to (tr(P 7 )) <~ 0.2 
for r g = 1.5. 

4. RESULTS AND DISCUSSIONS 

The convergence field, k, is estimated using the original 
Kaiser & Squires (1993) inversion algorithm with a Gaussian 
window function W(9) = exp(-|0| 2 / 0^) / 'ttOq. In order to esti- 
mate the noise of n field, we randomized the orientations of the 
galaxies in the catalog and created a K n oise map. The variance 
of the K no i S e depends on the smoothing scale G and is repre- 
sented as cr no i S e ~ O.O230 C '. This turned out to be consistent 
with theoretical estimate, 

- Ut <x\ 
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where a e = 0.42 and n g = 39.3arcmin~ 2 are the observed vari- 
ance of galaxy ellipticities and number density, respectively, 
yielding cr noise = O.O2O0 G '. 

We repeated the randomization 100 times and compute the 
rms among them at each grid where k is computed. Since k 
is Gaussian in its errors, this rms represents l-a noise level on 
each grid, and thus the measured signal divided by the rms gives 
the S/N ratio, v, at that grid. Thick contour in Fig. 1 shows the 
S/N-map with a smoothing scale of G = 1'. Threshold of the 
contour is set at 3. 

One notes that numbers of dark halo candidates are detected. 
For reference, number density of moderately bright galaxies are 
super-imposed as thin contour. It is surprising that correlation 
of mass and light is clearly visible from this relatively simple 
analysis on the single band image. A halo that has the high- 
est significance is designated (X) in Fig. 1. In this halo, the 
galaxy density also shows an excess as large as the other known 
clusters. In fact, a spectroscopic observation of two member 
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galaxies implies that it is a cluster located at the redshift of 
0.41. This object is not listed on the NASA/IPAC Extragalactic 
Database (NED). Other follow-up observations are underway 
and detailed comparison of mass and light will be discussed in 
subsequent papers. 
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FIG. 1. — Thick contour represents the signal to noise ratio (v) of 
convergence field. Data of v > 3 is shown. Thin contour shows num- 
ber density of moderately bright galaxies, 20 < R c < 23. Symbols (a)- 
(h) are designated where clusters listed in the NASA/IPAC Extragalactic 
Database (NED) are located and excess of galaxy number density is visi- 
ble in the map; (a)GHO 1606+4346, (b)Abell2158 (0.13), (c)GHO1601+4259 
(0.54), (d)GHO1601+4253 (0.54), (e)GHO1604+4303, (f)GHO1603+4256, 
(g)GHol602+4245, (h)GH01559+4242. The redshift is shown in the bracket 
if known. Both maps are smoothed with Gaussian kernel of 9 G = 1'. 

Let us turn to statistical properties of the S/N-map, and we 
will compare the results with the model predictions. First, we 
examine the probability distribution function (PDF) of peaks. 
Peaks are identified by pixels having a higher or lower value of 
v than the all eight surrounding pixels as suggested by Jain & 
van Waerbeke (2000). The PDF of the peak heights is shown as 
thick histograms of left two panels in Fig. 2. The thin histogram 
shows a noise PDF measured from K no i S e maps, the average of 
100 randomizations is plotted and the error bars show the rms 
among them. The bimodal feature of the noise PDF is due to 
superposition of two contributions each of which is due to pos- 
itive peaks and negative peaks (= troughs). Excess counts over 
the noise PDF is apparent on both positive (v > 3) and negative 
(y < -3) side of the histogram. The excess becomes more ap- 
parent as 9 G increases. The positive peaks imply the existence 
of dark matter halos whereas the natural explanation of the neg- 
ative peaks seen is voids (Jain & van Waerbeke 2000). This is 
the first observational evidence of the existence of voids in the 
matter distribution revealed by weak lensing. 

Right panels of Fig. 2 show the peak PDF measured from 
numerical experiments of the gravitational lensing by large- 
scale structures. The experiments were performed using the 
ray-tracing technique combined with large N-body simulations 
and its details are described in Menard et al. (2002). Briefly, 
Af-body data from Very Large simulation carried out by the 
Virgo Consortium (Jenkins et al 2001; Yoshida, Sheth & Di- 
aferio 2001) were used to generate the dark matter distribution. 
A ACDM model (O m = 0.3, f2 A = 0.7, h = 0.7 and cr 8 = 0.9) is as- 
sumed. The multiple-lens plane ray-tracing algorithm was used 
to follow the light rays (Hamana & Mellier 2001, see also Jain, 
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Seljak & White 2000 for technical detail). The lensing con- 
vergence and shear were computed for 1024 2 pixels with the 
pixel size of 0.25 arcmin for a single source plane of z s = 1. We 
generated the simulated noisy k map following the procedure 
described in Jain & van Waerbeke (2000). The noise is mod- 
eled as a Gaussian random field with the variance measured in 
the real data. The peak is identified as the same manner as for 
the observed data. Ten realizations of the numerical experi- 
ment were used to compute the peak PDF (the total field size is 
10 x 4.27 2 deg 2 , much larger than the observed data, and thus 
the peak distribution is very smooth). 
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FIG. 2. — Peak distribution function in kappa map with smoothing scale 
9 G = 0'.5 (top) and 1' (bottom). Thick solid histograms on the left panels 
show the results obtained from the observed catalog and thin solid histograms 
associated with error bars are averaged results of 100 random realizations of 
galaxy orientation in the catalog. The errors are rms of the randomization. 
Results from ACDM model simulation are shown as thick histograms on the 
right based on the observed parameters, (J noise 

= O.O230 G assuming z s = 1.0 
and their pure noise PDF is shown as thin histograms. 

The simulated PDF is very similar to the observed PDF. Jain 
& van Waerbeke (2000) pointed out that the peak PDF can be 
used to discriminate the cosmological model. They demon- 
strated using the ray-tracing experiments of a high density 
(n m = 1 .0, n x = 0, cr 8 = 0.9) and a low density(fi m = 0.3, fi A = 0, 
a& = 0.9) models that the high density model has broader dis- 
tribution than the low density model, especially on the negative 
side. It should be noticed that they adopted a e = 0.2 which is 
about half of our measured value, therefore one should not di- 
rectly compare our real PDFs with theirs (Figure 3 of Jain & 
van Waerbeke (2000)). 

In addition to the limited statistics of this work, theoretical 
prescription of the peak PDF is, however, still premature for 
placing constraint on the cosmological model in a quantitative 
manner. We therefore concentrate only on excess on the posi- 
tive side here; the halo number count. We calculated number of 
the positive peaks, N b s , in the observed S/N map exceeding z/ t h. 

Halo number count is defined as Nhaio = N h s -N r ± N b s + a\ ir 
where N r and <r Nr is average peak counts obtained on the ran- 
domized catalogs and their variances, respectively. The halo 
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count is shown in Fig. 3 as well as theoretical predictions which 
are constructed following Kruse & Schneider (1999). We as- 
sumed a high positive peak in k map comes from a single halo 
whose number density is described by Press-Schechter mass 
function (Press & Schechter 1974; we used its modified model 
by Sheth & Tormen 1999). We further assumed that all ha- 
los have a universal density profile that Navarro et al. (1996) 
proposed (NFW); p(r) = p s r~ x (\ + r)~ 2 . We also compared the 
results with another halo profile, singular isothermal sphere 
model (SIS); p(r) = p s r~ 2 . For the relation between the scale 
length of the halo and halo mass we adopted the fitting form by 
Bullock et al (2001). The cr noise measured from observed data is 
used to convert k to v. Note that Kruse & Schneider (1999) and 
Bartelmann et al. (2001) adopted cr e = 0.2, thus their predicted 
count is much larger than our prediction. 
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First, if one considers only highly significant halos, say v > 
5, the predicted halo number count does not depends strongly 
on the background cosmology as is presented in Fig. 3, where 
dashed line and solid line shows the SCDM model (Q m = 1, 
fl A = 0, cr 8 = 0.6, h = 0.5) and ACDM (Q m = 0.3, n A = 0.7, 
a% = 0.9, h = 0.7), respectively. The count, however, depends 
on the halo profile; the SIS predicts less counts as is shown in 
the dotted line in the figure. The observed count prefers NFW 
profile although the poisson error in v > 6 is large in this work. 
It is demonstrated here that number count of highly significant 
halos can be an useful probe for the dark matter halo profile 
because of the independence of background cosmology. Mean- 
while, the dependence becomes clearly visible in the counts of 
moderately significant haloes (v < 5); the SCDM predicts more 
counts than LCDM does due mainly to larger fl m . Our observed 
result shows high number counts which favors high density uni- 
verse. This result, however, cannot be taken at face value since 
the large excess is apparent in Fig .2 as a spike around ~ 3.5. 
Therefore, the excess count in the moderately significant ha- 
los is likely due to the cosmic variance. Comparison of data 
taken on multiple (say 5 to 10) uncorrected fields of this size is 
necessary to evaluate the contribution from the cosmic variance 
and eventually to discuss cosmological models. 

In summary, the present work clearly demonstrates that weak 
lensing technique is indeed a promising way to probe matter 
distribution purely through the mass. We detected not only the 
high positive peaks due to massive halos but also the excess in 
the peak PDF on the negative side which is the first observa- 
tional evidence of the voids in dark matter distribution. Dark 
matter halo counts is proved to be the useful tool to test the 
paradigm of the structure formation and cosmological models. 
The cosmic variance that we encountered in this work should be 
evaluated by the data obtained by future high resolution wider 
field imaging survey. 



FIG. 3. — Dark matter halo counts on the convergence field with 9a = 1'. 
The counts due to noise are estimated from the random realizations and sub- 
tracted. Solid and dotted line show ACDM model with the halo profile of 
NFW and SIS, respectively. Dashed line presents SCDM with NFW profile. 
The model calculation assumes a single source redshift of z s = 1 .0. Broadening 
of the redshift distribution makes the count slightly less but the effect is a few 
percents at most as far as the empirical distribution with reasonable parameters 
are assumed (e.g. Hamana et al. 2002). 
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